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The Inflationary Trispectrum for Models with Large Non-Gaussianities 
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We compute the leading order contribution to the four-point function of the primordial curvature 
perturbation in a class of single field models where the inflationary Lagrangian is a general function 
of the infiaton and its first derivative. This class of models includes string motivated inflationary 
models such as DBI inflation. We find that the trispectrum for some range of parameters could 
potentially be observed in future experiments. Moreover, the trispectrum can distinguish DBI 
inflation from other inflation models with large non-Gaussianities which typically have a similar 
bispectrum. We also derive a set of consistency conditions for n-point functions of the primordial 
curvature perturbation in single field inflation, generalizing Maldacena's result for 3-point functions. 



The Cosmic Microwave Background (CMB) provides 
us with a remarkably detailed snapshot of the early uni- 
verse. The vast structure that we see today is consistent 
at present with the simple picture of an almost scale in- 
variant, Gaussian primordial density perturbation gener- 
ated by inflation. Precision measurements of any small 
deviation from this simple picture contain valuable infor- 
mation that could constrain cosmological models. 

In the case of single field inflation, the non-Gaussian 
perturbation is governed by the n-point functions (n > 
3) of the curvature perturbation £. The leading non- 
Gaussian features are known as the bispectrum (three- 
point function) and trispectrum (four-point function), 
with their sizes conventionally denoted as /jvl, t nl re- 
spectively. The current experimental bound for the bis- 
pectrum from WMAP3 is -54 < f NL < 114 @. Al- 
though the experimental bound for the trispectrum is 
rather weak at the moment |tjvl| < 10 s B, the next 
generation of experiments such as PLANCK will increase 
the sensitivity to about \tnl\ ~ 560 j|, see also [fl6"| , 
pp| . The CMB bispectrum and trispectrum, if observed, 
could provide strong constraints on inflationary models. 
In particular, the shape of the trispectrum contains even 
more information about the inflationary dynamics than 
the bispectrum because of the richer varieties of momen- 
tum dependence it can encode. As we shall see, in some 
string motivated inflationary models such as DBI infla- 
tion pj , the size of the bispectrum and trispectrum can 
both in principle reach the observable limit of future ex- 
periments. Therefore it is crucial to study the shapes of 
higher-point functions in such inflationary models. 

The CMB bispectrum for single field inflation has been 
studied intensively, see, e.g., (§, |L @, The four- 

point function for slow roll inflation, however, has only 
recently been studied [jl(|, where it was found that the 
inflationary trispectrum generated by slow roll models 
is too small to be observed even by future experiments. 
In this note we compute the shape of the trispectrum for 
DBI inflationary models and more generally a class of sin- 
gle field models considered in jllj] where the inflationary 
Lagrangian is a general function of the infiaton and its 
first derivative. We find that for n-point functions, there 



are n — 1 independent shapes that could be potentially 
large, providing useful information in distinguishing vari- 
ous inflation models. In the context of slow-roll inflation, 
Maldacena || has proposed an interesting consistency 
relation for the three point function in a squeezed limit 
where one of the momentum modes has a much longer 
wavelength than the other modes. The importance of 
this relation is that its violation, if observed experimen- 
tally, could rule out a large class of inflationary scenarios 
model-independently |l7| . In view of this, we derive a set 
of consistency conditions for higher-point functions, gen- 
eralizing Maldacena's result. Although some aspects of 
the consistency conditions for the trispectrum in slow roll 
models have recently been discussed in jH^, our result 
is applicable to a more general inflationary Lagrangian 
and to cases when there are more than one long wave- 
length modes. We show that the shapes of the trispec- 
trum we found satisfy these generalized consistency con- 
ditions. We end with some final remarks and discuss the 
implications for future CMB experiments. 

Let us begin with DBI inflation which was proposed 
in [Q and studied subsequently in e.g. @, |l2|, |l|, 
Although DBI inflation can be formulated without re- 
ferring to branes and extra dimensions, it is most easily 
visualized as describing a D-brane moving relativistically 
in a warped throat where the infiaton <fi is the position 
of the D-brane. For a warp factor /(</>) and an infiaton 
potential V (((>), the inflationary Lagrangian is 



S = 
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is the Hubble 



where a(t) is the scale factor and H = ^& 
parameter in the Friedmann-Robertson- Walker (FRW) 
universe. The background solution for the infiaton <f> 
is spatially homogeneous. The sound speed, which de- 
scribes the speed of fluctuations on this background, is 
defined as 



i - f{4>W 



(2) 



If the potential V(tf>) is not fiat, the inflaton will quickly 
approach the relativistic speed limit cb 2 ~ /(</>) _1 and so 
c s <C 1 • We will see that for DBI inflation, the bispec- 
trum and trispectrum scale with the sound speed as: 
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4 denote the orders of a in these functions. We use the 
definition of the momentum density 



dC 

da 



(6) 



(3) to express a in terms of it up to the third order, 



This can be easily understood as follows. By expanding 
the DBI Lagrangian in the small sound speed limit, we 
find that the n-point function scales like . Since /jvl 
is proportional to the three-point function divided by the 
square of the power spectrum, it scales as ^. Similarly, 
t nl is proportional to the four-point function divided by 
the cube of the power spectrum, and hence it scales as 
-TT. Thus, there can exist models with a sufficiently small 
sound speed c s that the size of the trispectrum tnl is 
bigger than the sensitivity level of future experiments 
t nl > 560, but yet not too small to be already ruled 
out by the WMAP bound on the CMB bispectrum f NL . 
As analogous to the bispectrum, we expect the bound on 
tnl to depend on its shape [[t5| and so it would be inter- 
esting to carry out a similar analysis for the trispectrum. 
We note, however, that the bound on the trispectrum is 
rather weak and so only in the small sound speed limit 
and the leading contributions to tnl could possibly be 
observable in future experiments. In the small sound 
speed limit, the non-Gaussianities are dominated by the 
inflaton's derivative interactions. The contributions from 
the gravity sector and potential terms of the scalar are 
suppressed by slow roll parameters, we will not consider 
these subleading contributions in this paper. However, 
it is of theoretical interests [|l(| to study in more details 
these subleading contributions since they may provide a 
testing ground for theoretical issues such as the gener- 
alized Maldacena's consistency conditions that we will 
present here. 

We use the "in-in" formalism jll| to calculate the four- 
point expectation value for the primordial curvature per- 
turbation, 

<C 4 ) = -* f dt / ([C I (t,k 1 )--< I (t,k 4 ),H int (t')}) , (4) 



where H mt — J d 3 xrii n t is the interaction Hamiltonian. 
Consider a perturbation a(t, x) = 6cj> around the back- 
ground solution. We first work out the general relation 
between the interaction Hamiltonian and the interaction 
Lagrangian up to the fourth order in a. We expand the 
fluctuation Lagrangian density (starting from the second 
order) to the fourth order in a(i,x), and denote it as 

£ = foa 2 + h + .9o« 3 + 9ia 2 + 92& + h 

+ hpa 4 + /iid 3 + h 2 a 2 + h 3 a + ji , (5) 

where the fp, j's, g's and h's are all functions of a(t, x), 
its spatial derivative did and t. The subscripts 0, 1, 2, 3, 
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and plug it into the Hamiltonian density 
7i = ira — C . 



(7) 
(8) 

(9) 
(10) 



We then separate the TL into a kinematic Hamiltonian 
density Hp which is quadratic in -k and a, and an in- 
teraction Hamiltonian density 7ii n t- To use the "in-in" 
formalism, the a and 7r in this interaction Hamiltonian 
should be replaced, respectively, by the a 1 and it 1 in the 
interaction picture. These latter variables satisfy the free 
equation of motion (in the time-dependent background) 
followed from Tip and the usual commutation relation 



[a I (t,x),TT I (t,y)]=i6 3 (x-y) . 
Expressing the ir 1 in 7i j n t in terms of a 1 using 
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we finally get (omitting the label "I" in the variables in 
Hi n t from now on) 
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As we see, while the cubic part of 7ij n t is the negative of 
that of £i n t , this is generally not so for higher order ex- 
pansion. For slow-roll inflation, the extra terms are sup- 
pressed by slow-roll parameters |Q, and therefore can 
be properly neglected in the trispectrum calculation of 
Rcf. |l(]]. As we will see, for the case of interest here, the 
extra terms contribute to the leading order results. 

For the special case of the DBI action, to compute the 
leading shape of the four-point function, we only keep 
the leading terms in the small sound speed limit, which 
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come from expanding the square root term in the DBI 
action. The Lagrangian density up to quartic order is: 



£2 = —^-a^ci^JaY] , 
£ 3 = _^! T [ £ j3„ a -2 c 2 d(Va) 2 
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Using (113) we then have 
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(15) 
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From the quadratic term we can solve for the perturba- 
tion and quantize it according to the standard procedures 
of quantum field theory: 



a(t, x) 



1 



d 3 k[u(r, k)a(k) 



(2tt) 3 

+u*(T,-k)o t (-k)]e ik - x , (19) 

where u(t, k) = -|= ( 1 + ikc s T)e~ lkCsT is the solution of 

the quadratic Lagrangian, and t = — is the conformal 
time. The creation and annihilation operators satisfy the 
usual commutation relation [a(k), a^(k')] = (27r) 3 <5 3 (k — 
k'), consistent with (pi]). 

In the small sound speed limit, the primordial density 
perturbation £ is dominated by the fluctuation of the 
inflaton a, and is simply related to it by C = ^tO- @- 
The two point function is 



<C(k 1 )C(k 2 )) = (27r) 5 ( 5 3 (k 1 +k 2 )^, 



(20) 



where P k = 



is the power spectrum, and we have 
evaluated the perturbation at several e-foldings after the 
horizon exit with the approximation r — > 0. 

Using (|l8|), the four-point function at leading order in 
C is 

(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) 

= (2^(^k i )(^) 3 ^ n ^3(^ 1 ), (21) 



where 



Al = -72HiM 
1 K 5 



(22) 



Here P^- is the power spectrum at the leading order in 
slow variation parameters. Momentum conservation im- 
plies that the momenta ki where i = 1 . . . 4 form a quadri- 
lateral. The quadrilateral may be skew, i.e. not lying 



on a plane. For convenience we will just refer to it as 
a quadrilateral. This is in contrast to the bispectrum 
where the momenta form a triangle. 

In order to compare with the value of tn l from exper- 
iments, we need to specify a configuration of the quadri- 
lateral to evaluate the four-point function. We consider 
the "equilateral configuration" in pJ, where the 4 mo- 
menta have the same amplitudes k, and hence the an- 
gles between them satisfy cos(#i 2 ) = cos(# 34 ), cos(#i 3 ) = 
cos(6» 24 ), cos(6»i 4 ) = cos(6> 23 ), and Y? m =i cos(0 m4 ) = -1. 
Here Oij is the angle between k^ and k, . Using the defi- 
nition of tnl 

(C(k 1 )c(k 2 )C(k 3 )C(k 4 )) 



tnl 



(2^) 9 5 3 £>) 



x [P c (A:i)P c (fc 2 )P c (fci 4 ) + 23 permutations] , (23) 
where Pc(fci) = P^/(2fe?) as defined in (|20|) and k, 



k, + k 

and 



mj - 

3 , we find that for DBI inflation, by matching (|2] 
|2"3|) in the special "equilateral configuration" , 



2§ 



tnl 



128 ^E; l= i(i+cos^ m4 ) 



(24) 



For a generic angular configuration, we can roughly esti- 
mate that tnl ~ ^ji - - This confirms our claim that there 
are models that give rise to observable trispectrum, but 
are not ruled out by current bounds on the bispectrum, 
e.g., when c s ~ 0.05. 

It is straightforward to generalize the analysis to a class 
of models where the inflation Lagrangian is a general 
function P(X,<p), with X = -\g» y d^(j)d u cj> (ll) . The 
sound speed is defined as 



<?.= 



P 



.X 



Px + 2XP 



XX 



(25) 



Again we consider only the leading contribution in the 
small sound speed limit. In this case the derivative of 
the Lagrangian with <j>, i.e., P^, does not give rise to 
shapes large enough to be observed. One can show that 
the leading order expansions up to quartic order are 



C l ( dP 
C2 2 { dX 
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(27) 
a~ 2 {Vaf) 
(28) 



So for the four-point function we found 3 independent 
shapes (i.e., momentum dependences). A convenient ba- 
sis is the shape A% we found earlier, and two additional 
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shapes coming from terms proportional to ad 2 (Va) 2 and 
a _1 (Va) 2 (Va) 2 , respectively, in Ti 
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.4, 



lfc 2 fc 2 (k 3 -k 4 ) 
V 
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3(k 3 + fc 4 ) I2fc 3 fc 4 
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K 2 

k 1 k 2 k 3 k 4 
K 1 ' 



K 2 



(29) 



(30) 



where K = J2i=i anc ^ "perm." refers to the other 23 
permutations of the four momenta. The coefficients of 
these shapes are determined by the Lagrangian (|26|)-(|2T 
and the relation (|l4|). 

The shapes of the four-point function A\, A 2 , A 3 can 
be used to distinguish between models with different pa- 
rameters §jf[, §jfs, §7T*- ^o illustrate this, we can com- 
pare the situation with that of the bispectrum, where 
the shapes were computed in fA 61. There it was found 
that at the leading order there are two large shapes of 
non-Gaussianities, denoted by A\ and A c : 



A x 

A c . 



k 2 k 2 k 2 

K 3 ' 

i>j 



rr k- k ■ 

2K 2 ^ 1 3 



X 3 -(31) 



However, these two shapes are functions of the momen- 
tum triangle and look qualitatively similar in plots . In 
particular, in the squeezed limit k\ — > 0, one can easily 
show these two shapes have the same asymptotic behav- 
iors A\ , A c ~ kj. The shapes of the four-point function 
Eqs. j22]) , (p9|) , (|30|) are more distinguishable. For exam- 
ple, we can take a squeezed limit k\ — > 0, and the re- 
maining three sides of the quadrilateral form a triangle. 
The asymptotic behaviors of the shapes are 
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[/c 3 fc 4 cos(6»i 2 )(l 
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/c 2 fc 2 + fc 2 fc 4 + &3&4 3k 2 k 3 k 4 



K v ' ^2 A 3 ' 

[cos(0i 2 ) cos(6» 34 ) + cos(6»i 3 ) COS(0 2 4) 

+ cos(0 14 )cos(0 23 )] . (32) 



In the squeezed limit of a generic configuration of the 
quadrilateral, the first shape scales as A 4 ~ A: 2 , and two 
other shapes scale as A 2l A 3 ~ fci. So this in principle 



provides a more refined discriminator between different 
inflation models than the bispectrum. Another interest- 
ing configuration is when ki, k 2 and k 3 are orthogonal 
to each other and k 4 = — ki — k 2 — k 3 . In this configu- 
ration, the shape A 3 vanishes. In particular, these pro- 
vide the distinction between DBI inflation and the power 
law k-inflation model in The simplest power law 

k-inflation has a Lagrangian P(X, <fi) ~ (—A + X 2 )/4> 2 . 
It was shown in jOJ |l) that the model has a power law 
inflationary solution and the inflaton moves at constant 
speed. Using the interaction Hamiltonian (|l4|), we find 
the fourth order Hamiltonian has three terms propor- 
tional to a 



3 <9vP • 4 



f.V- 



a , a 



P A,2 



a (Va) , a 



-i|^ (V a) 2 (Va) 



We know for example from the calculations in the DBI 
case that the d 4 , a~ 2 c 2 a 2 (ya) 2 7 a~ 4 c 4 (Vo;) 2 (Va) 2 give 
rise to the same c s order contribution to trispectrum 
in the small sound speed c s limit. So we can see the 
A 3 shape coming from the interaction term (Va) 2 (Va) 2 
dominate over the other shapes for the simple power law 
k-inflation. On the other hand, the DBI inflation does 
not have the A 2 and A 3 component and their features 
in the above momentum configurations. So these distin- 
guish the power law k-inflation from the DBI inflation. 

One can go to more complicated corners of the config- 
uration space of the quadrilateral and study the different 
behaviors of the shapes A±, A 2 , A 3 , but it should be 
clear that the configuration space of a quadrilateral is 
much bigger than that of a triangle, and thus provides a 
more powerful discriminator for inflationary models. 

Now we derive a generalized consistency condition for 
rt-point function, and test it in various models for which 
we have computed the four-point function. In single 
field inflation, the density perturbation is generated by a 
scalar perturbation £ whose 1-point function is: 

(C(ki) ■ ■ ■ C(k ; )> = (27r) 3 ,5 3 (k 1 + ■ ■ ■ + ki)P ! (ki, ■ ■ ■ , kj) (33) 

The scalar perturbation ( has scaling behavior of degree 
—3, so the /-point function P l (ki , • • • , kj) is roughly a de- 
gree 3 — 3Z homogeneous function of ki, i = 1, • • • , I. Wc 
can define the spectral index of /-point functions analo- 
gously to that of the two-point function 



ni 



I 



dlog(P'(ki 



<flog(fc) 
eflog^ki)---^))) 



3/ - 3 



•31 



(34) 



d\og{k) 

In the above definition all the kj are taken to be of the 
same order of magnitude as k, e.g. = 1. When 

/ = 2, this is just the usual spectral index of the power 
spectrum n s — 1. Another example is 1 = 3. In this case 
the three-point spectral index is 

d\og(f NL ) 



n 3 -l = 



dlog(fc) 



+ 2(n s - I) 



(35) 



where dlo s(fNL) 
Gaussianity p2j 



is known as the running of non- 
In an inflationary model where the 
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physical parameters vary slowly within a Hubble time, 
the /-point spectral index is very small, i.e. \ni — 1| <C 1- 
We will take the squeezed limit of the n-point function, 
where I of the momenta (2 < I < n) are much bigger than 
the others, i.e., 



as , while the right hand side scales as j^-ts 



so 



fell 



, h > ki+i, ■■■ ,k n 



(36) 



The long wavelength modes fcj+i, • • • , k n exit the horizon 
and are frozen much earlier than the short wave modes 
ki , • • • , ki . Their effect is to act as a background that 
scales the spatial coordinates, since they perturb the 
FRW metric as ds 2 = dt 2 - a{t) 2 e 2 ^ B dx 2 . The /-point 
function in a background perturbation C,b ~ 10~ 5 is 

(C( Xl ) • • • C(x ; ))b = (C(e CB xx) • • • C(e Cs x ; )) . (37) 

The background Cb is roughly constant around the short 
wave length scale \xi — Xj\ with 1 < i < j < I, and can 
be evaluated at x = Xl +y+ x ' . We expand the above 
equation to first order in the background 

(C(xi)---C(x,)) B = (C(xi)---C(x I )> 
+CB(x ){(^-(C(e Cfl Xl ) • • • C(e Cfl Xi )))| CB =o} (38) 

The zeroth order term is independent of the background 
so will not contribute to the final result. After a Fourier 
transformation to momentum space, the first order term 
in the background expansion becomes 



d 3 kl 



CB(kB){ 3 ^(e- 3/CB (C(e- Cfl (k 1 



(2tt)3^— d( B 
.. C (e-^(k,-^))>)| te=0 } 



-)) 



(39) 



Here we introduce an integration variable k^, which will 
be fixed by the delta function to be k^ = ki + • •• + k/, 
and which satisfies ks <C ki, i = 1, • • • I. At the leading 
order of the squeezed limit we can neglect the background 
momentum in the /-point function, and use our definition 
of the /-point spectral index to find 

(C(ki)-(W}B = -(^-^(kr,---,^) 

xC(k 1 + --- + k / ) (40) 

Plugging in the expression for the n-point function, we 
find a generalization of Maldacena's consistency condi- 
tion to n-point function in the squeezed limit k\ , ■ ■ ■ , fe; S> 
ki+i, ■ ■ ■ ,k n : 



P n (k 



1) • • " 



,k„) = -(n / -l)P / (k 1 ,---,k / ) 



xP"- ,+1 (k 1 + ...k i ,k w ,...,k„) (41) 

In the case of DBI inflation, the n-point function on 
the left hand side of the consistency condition (fhij) scales 



the n-point function must vanish in the squeezed limits 
at leading order for small sound speed. We see indeed the 
shapes A 1: A 2 and A 3 in (||), (||), (|o|) all vanish in a 
squeezed limit when there is a long wave mode ki — » 0, in 
agreement with the n-point consistency condition ( |4l"| ) . A 
more non-trivial check of (41 ) is when P n (ki, • • • , k n ) is 
non- vanishing in the squeezed limit, e.g., one may check 
the validity of (^) using the recent results for slow roll 
inflation jl^] (together with subleading terms which are 
yet to be determined) . Besides serving as a check on the 
calculations, the generalized consistency condition ( (jd] ) if 
violated can, in the spirit of |l7|| , rule out a large class of 
models experimentally. 

For a general Lagrangian P(X,(f>), the n-point func- 
tion has n — 1 independent shapes proportional to the 
derivatives of the Lagrangian J^r, where i = 2, 3, • • • ,n. 
It is not difficult to compute the shapes and show that 
they vanish in the squeezed limit, though it is difficult to 
measure such higher point functions experimentally. 

We emphasize that the CMB trispectrum, although 
difficult to observe, could provide a useful discrimina- 
tor for inflationary models. In particular, if large non- 
Gaussianities are observed experimentally in the future, 
it may still be difficult to distinguish via the bispectrum 
DBI inflation from other inflationary models governed 
by derivative interactions, e.g., the models in J9|. The 
trispectrum we computed may serve as an extra filter to 
zero-in on string inflationary scenarios. 
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